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Comment on "Polynomial-Time Simulation of 
Pairing Models on a Quantum Computer" 

In a recent Letter Wu et al. proposed a polynomial- 
time algorithm (PTA) for simulations of a wide class 
of pairing hamiltonians on a NMR quantum computer. 
While this is an interesting result we want to recall in 
this Comment the existence of classical algorithms for 
attacking this problem to any required precision. There 
are two main sort of algorithms. If the pairing model 
(PM) turns out to be integrable then it can be solved a la 
Bethe. This happens in the reduced BCS model charac- 
terized by a single coupling constant, but there are more 
general cases reviewed below. The other algorithm is the 
density matrix renormalization group (DMRG), which is 
applicable to both integrable and nonintegrable models. 

The basic idea of the Letter is that within the next 
generation of quantum computers (QC) with up to 10-100 
qubits, it will be possible to solve hard quantum many 
body problems in Hilbert spaces with N ~ 10 2 single 
particle levels and M ~ 10 2 number of particles. However 
the authors seem to have overlooked the aforementioned 
facts concerning the power of analytical and numerical 
methods as the DMRG. While the analytical methods 
are only applicable to integrable models, the DMRG is 
of general applicability. 

As mentioned in [l| the PM is of fundamental impor- 
tance for nuclear structure studies in medium and heavy 
nuclei and it was introduced in nuclear physics 0] soon 
after the work of Bardeen, Cooper and Schieffer. Even 
though the BCS theory explains qualitatively the super- 
conducting phenomenon in finite nuclei, it suffers from 
the strong number fluctuations. 

More recently the (PM) Hamiltonian has been exten- 
sively used to describe the physics of ultrasmall supercon- 
ducting grains. It was shown in this context that only by 
resorting to the exact numerical solution 0], using the 
DMRG, was it possible to describe the disappearance of 
superconducting correlations as a function of the grain 
size. The ground state energy for systems as large as 
N ~ 400 with Hilbert space dimensions of ~ 10 130 were 
treated with high numerical accuracy (5 significant fig- 
ures). In particular, the case of N = 100 was selected to 
study the convergence properties of the method showing 
that the exact ground state energy can be obtained to 
seven significant figures with a moderate computational 
effort (see Table I of Ref. @|). The DMRG is a highly 
accurate iterative method based on the Wilson renormal- 
ization group. In the infinite algorithm procedure used 
in Q the number of iterations needed to obtain the low 
lying states of a pairing Hamiltonian is N/2 which has 
to be compared with the ~ N computational steps re- 
quired by PTA of Ref. ij. At each iteration de DMRG 
procedure stores the matrix elements of all relevant oper- 
ators in an optimum subspace of m selected states. Due 
to the exponential convergence of the method m is tuned 


to obtain the desired precision. For a general pairing 
hamiltonian there are three operators per level, therefore 
the total memory usage is (3(3m 2 N) which for reasonable 
values of m ~ 10 2 allows the treatment of systems with 
N ~ 10 3 . In we used a constant pairing interaction 
but the DMRG approach can easily accommodate arbi- 
trary pairing matrix elements, and it is specially suited 
to give the ground state and the low lying excited states. 

Immediately after this work appeared, it was realized 
that the PM has been solved exactly by Richardson in the 
1960's 0. Recently 0, the exact solution of the PM has 
been generalize to a large class of pairing Hamiltonians. 
The most general pairing Hamiltonian can be written as 


where 


- c f 


and i runs over the N 


double degenerate levels. Note that the third term 
corresponds to a monopole interaction absent in the 
Hamiltonian treated in the Letter. As shown in ref. 

there is a wide subset of exactly solvable Hamil- 
tonians parametrizing the interaction as £j = u, — 


Efcf^il 7 (e; - <y) cot 7 ( Vl - rjj), V,} 


2g 


7(ei~gj) 


and V? = §7 (e, - Cj) cot 7 (77, 


Cj) cot 7 (rji — T)j). Where g, and 
r/i are a set of 2N + 1 free parameters and 7 can take the 
values 0, 1 and — i defining the three families of exactly 
solvable models, the rational, the trigonometric, and the 
hyperbolic models respectively. As we see, among the 
2A 2 — N set of free parameters of the most general and 
non-integrable Hamiltonian i ffl l, a restriction to a sub- 
set of 6N + 3 free parameters leads to an exactly solv- 
able Hamiltonian. We believe that most of the physi- 
cal problems can be modelled with a pairing Hamilto- 
nian within the integrable subset, but if some particular 
case requires a non-integrable pairing Hamiltonian, it can 
then be solved numerically to the desired precision using 
the DMRG. 
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